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Abstract 

In the present work, we provide a dielectric study on two differently concentrated aqueous lysozyme solutions in the 
frequency range from 1 MHz to 40 GHz and for temperatures from 275 to 330 K. We analyze the three dispersion 
regions, commonly found in protein solutions, usually termed /3-, y-, and ^-relaxation. The /^-relaxation, occurring 
in the frequency range around 10 MHz and the y-relaxation around 20 GHz (at room temperature) can be attributed 
to the rotation of the polar protein molecules in their aqueous medium and the reorientational motion of the free 
water molecules, respectively. The nature of the 5-relaxation, which often is ascribed to the motion of bound water 
molecules, is not yet fully understood. Here we provide data on the temperature dependence of the relaxation times 
and relaxation strengths of all three detected processes and on the dc conductivity arising from ionic charge transport. 
The temperature dependences of the fi- and y-relaxations are closely correlated. We found a significant temperature 
dependence of the dipole moment of the protein, indicating conformational changes. Moreover we find a breakdown 
of the Debye-Stokes-Einstein relation in this protein solution, i.e., the dc conductivity is not completely governed by 
the mobility of the solvent molecules. Instead it seems that the dc conductivity is closely connected to the hydration 
shell dynamics. 

Keywords: protein dynamics, protein solutions, dielectric spectroscopy, lysozyme, hydration shell, 
Debye-Stokes-Einstein relation 



1. Introduction 

Proteins are essential for life. They are the build- 
ing blocks of cells and they are part of virtually ev- 
ery biological process IH As enzymes they cat- 
alyze chemical reactions and in cell membranes they 
build ion channels and pumps; they are responsible 
for signal generation and transmission and they also 
act as antibodies, hormones, toxins, anti-freezer, elas- 
tic fibers or source of luminescence. This ubiquity 
of proteins since long has triggered scientists' demand 
for a deeper understanding of their structure and func- 
tionality. In organisms, proteins with biological func- 
tions usually exist in solution and many of their physi- 
cal and functional properties are strongly influenced by 
the solvent |3]. Therefore it is vital to examine pro- 
teins within their common environment. Due to its im- 
portance for obtaining a deeper understanding of bi- 
ological processes, the dynamics of proteins ingen- 
eral is a very active field of research 10. S S 



which is often focused on the protein-water interaction 
12111311 . Here a suitable and commonly em- 
loyed experimental method is dielectric spectroscopy 
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Dielectric spectra of aqueous protein solutions show 
at least three dispersion regions y , 25 , H 27 ], revealing 
the typical signatures of relaxation processes, namely a 
step in the dielectric constant, e'(v), and a peak in the 
dielectric loss, e"(v) (Fig. [T]i. In the biophysics commu- 
nity, they are often termed y-, and 6- relaxation. This 
nomenclature should not be confused with that used in 
the investigation of glassy matter, where the terms a-, 
y6-relaxation, etc. are commonly applied to completely 
different phenomena than those considered here (see, 
e.g., refs. 28 and|H). In the present work we follow 
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the biophysical nomenclature. 

The y6-relaxation in the low frequency range and 
the y-relaxation at around 18 GHz (at room tempera- 
ture) can unambiguously be assigned to the rotation of 
the polar protein molecule in its aqueous medium and 
the reorientational motion of the free water molecules 
(similar to the main relaxation process in pure water). 
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Figure 1 : Schematic view of the dielectric loss spectram of a typical 
protein solution close to room temperature (circles). The shaded areas 
show the contributions of the /}, S, and y relaxations, which arise from 
reorientational motions of the protein molecules and the bound and 
free water molecules, respectively. The dotted line indicates typical 
raw data, dominated by dc conductivity at low frequencies. 



respectively 27]. (The y-relaxation corresponds 
to the ff-relaxation of bulk water within the glass- 
physics nomenclature. Within this nomenclature, the 
jS-relaxation may be regarded as a-relaxation of the pro- 
tein molecules, governed by the solvent dynamics.) The 
third dispersion, located between /3- and y-relaxation 
is still a subject of discussion. After the pioneering 
works of Oncley revealed the presence of (3- and y- 
relaxations in the late 1930's and 1940's HI [l5l[30ll, 
first indications for this 5-dispersion where found by 
Haggis and Buchanan about one decade later l3lll32ll . 
It is quite generally accepted nowadays that this disper- 
sion, which was detected in differ ent p rotein solutions, 
is due to bound water relaxation iiJO Is l Isllssll . 
It is well known that proteins possess a hydration shell 
of bound water molecules and it is reasonable that 
these water molecules should have slower dynamics 
than free molecules. However, the complexity of pro- 
teins makes it difficult to decide if the 5-dispersion can 
solely ll36ll be explained by a bound water relaxation 
(which could also be bimodal iHIl H H) or if ad- 
ditional effects hke intra-protein motions have to be in- 
cluded 1 18 , 38 , 2 , 24 , 2^ . The bound water relaxation 
has also been discussed in the context of the glass tran- 
sition in proteins 1221 • 

At low frequencies (< 1 kHz), the spectra of aque- 
ous protein solutions are dominated by electrode polar- 
ization (EP) lfl4l l23i I40I1 giving rise to giant values of 
the dielectric constant and a strong drop of conductivity 
towards low frequencies. EP arises when the conduct- 



ing ions in the sample arrive at the metallic electrodes 
and accumulate in thin layers immediately beneath the 
sample surface forming a so-called space-charge region. 
However, this effect is not a special feature of protein 
solutions, but affects dielectric spectra of any material 
containing free ions 1 4^ 41 , 42 ] . For this reason, here 
we will not present the low frequency region of our 
broadband spectra, dominated by EP contributions (see 
ref. 114011 for a detailed treatment of EP, including a 
lysozyme solution) . 

In the present work, we provide a thorough dielec- 
tric characterization of the relaxational processes in 
lysozyme solutions. Lysozyme, a representative of 
globular water-soluble proteins, is an enzyme and part 
of the innate immune system with a molar weight of 



14.3 kDa 1431 144|1- The obtained spectra in the fre- 
quency range from 1 MHz - 40 GHz allow for the de- 
tection of jS-, y-, and (5-relaxation. For the first time, we 
investigate the temperature dependence of spectra cov- 
ering all these intrinsic relaxations. This allows gather- 
ing valuable information including, e.g., the hindering 
barriers for the involved molecular motions. Further pa- 
rameters as the dipolar moment and the radius of the 
protein are deduced from the dielectric results and the 
validity of the Debye-Stokes-Einstein formula in this 
protein solution is checked. 

2. Materials and methods 

The complex dielectric permittivity and conductiv- 
ity were determined using two different coaxial re- 



In the frequency range 



flection techniques 1145 , 
1 MHz - 3 GHz, an Agilent Impedance/Material An- 
alyzer E4991A, was employed. The ac voltage is ap- 
plied to a platinum parallel-plate capacitor containing 
the sample material (diameter 4.8 mm, plate distance 
0.1 - 0.85 mm). The capacitor is connected to the end 
of a coaxial line, thereby bridging inner and outer con- 
ductor Contributions of the coaxial line and connec- 
tors were corrected by a calibration with three stan- 
dard impedances. For temperature-dependent measure- 
ments the capacitor is mounted into a N2-gas cryostat 
(Novocontrol Quatro). The sample holder and coaxial 
line, which connects the sample within the cryostat to 
the measuring device, were designed in our laboratory 
iH. The high-frequency range (100 MHz - 40 GHz) 
was covered by the Agilent "Dielectric Probe Kit" us- 
ing an open-ended coaxial line, the so-called "Perfor- 
mance Probe", in combination with a Agilent E8363B 
Network Analyzer Here, the line is immersed into the 
sample liquid, which is kept in 50 ml plastic tubes. The 
temperature was controlled by means of an Eppendorf 
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"Thermomixer Comfort" in combination with a 50 ml 
"Thermoblock" that is mounted on the heating plate of 
the Thermomixer to heat the sample tube. 

Dialyzed and lyophilized Lysozyme powder from 
chicken egg white was purchased from Sigma-Aldrich 
(Fluka 62970) and used without further purification. 
Aqueous protein solutions were prepared by dissolv- 
ing weighed amounts of protein powder in deionized 
H2O (Merck "Ultrapur"). The investigated concentra- 
tions correspond to 3 mmol and 5 mmol of Lysozyme 
powder added to one liter of water (equal to 42.9 mg and 
71.5 mg per 1 ml of water). The pH values of these so- 
lutions are around 3.8 (measured with a pH tester from 
Hanna-Instruments) . 

3. Results and discussion 

lysozyme smmoi 
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Figure 2: (a) Dielectric constant, (b) dielectric loss (coiTected for the 
contiibution from dc conductivity), and (c) real part of the conductiv- 
ity as function of frequency, measured at different temperatures. The 
lines are fits using the sum of a Debye function for the yS-relaxation 
and two Cole-Cole functions for the S- and y-relaxations. Inset: mag- 
nified view of the dielectric loss (b) in the region of /3- and (5-relaxation 
for five temperatures. 

Figure|2]shows the spectra of dielectric constant s'(v) 
(a), dielectric loss e"(v) (b), and real part of conductiv- 



ity cr'(v) (c) of a 3 mmol lysozyme solution covering 
the frequency range from 1 MHz - 40 GHz for three se- 
lected temperatures. The inset shows a magnified view 
of the loss in the region of the y6- and ^-relaxations for 
five temperatures in a semilogarithmic plot. 

In the real part of the conductivity, a strong 
dc-contribution is found showing up as frequency- 
independent plateau from 1 kHz (not shown) to 
100 MHz (Fig. |2|c)). This is due to ionic charge 
transport with the ions arising from the residual salt 
content in the protein sample (mainly chloride ions, 
left over from production process). Dielectric loss (b) 
and conductivity (c) are directly correlated via cr'(v) - 
e"{v)(DS() (oj = Ittv is the circular frequency and sq is 
the permittivity of vacuum). Thus, the dc conductiv- 
ity gives rise to a contribution s'^^ = crjf./(weo), i.e. a 
1 / V divergence in the loss, which obscures the detection 
of possible relaxation processes at low frequencies (cf. 
dotted line in Fig. [T]i. Therefore it is common prac- 
tice to subtract the dc contribution Js, 22, 24, 3^ 471. 
This is a critical task as the amplitude of the resulting 
relaxation peaks can strongly depend on the value of 
the subtracted dc conductivity. In the present case, the 
correctness of the subtracted values, which were deter- 
mined from the measured conductivity, is confirmed by 
the high quality of simultaneous fits of the step in s'{v) 
(unaffected by the dc conductivity) and of the peak in 
e"(v) revealed after subtraction (see below). This leads 
to the corrected loss spectra shown in Fig. |2jb). In 
this figure, clear signatures for the /3- and y-relaxation 
are found. They show up as peaks close to 10 MHz 
or 20 GHz, respectively, shifting to lower frequencies 
with decreasing temperature. This temperature depen- 
dence directly mirrors the reduction of reorientational 
mobility of the protein and free water molecules when 
temperature is lowered. However, a closer inspection 
of Fig. |2|b) provides clear indications for a third re- 
laxation process (5-relaxation) in the frequency range 
around 100 MHz. There is significant excess intensity, 
not explainable by a simple superposition of /3- and y- 
peaks. At the highest temperature shown (313 K), there 
is even the indication of a separate weak peak (see also 
inset of Fig. |2]i. However, clearly the y-relaxation is the 
dominating process and thus the and ^-relaxations are 
hardly discernible in the real part of the dielectric con- 
stant (Fig. |2la)). For low frequencies, e'(v) approaches 
a plateau whose absolute value of about 80-90 is of the 
same order of magnitude as the static dielectric con- 
stant of pure water (b' = 80.3 at 293 K H). The y- 
relaxation is sufficiently strong to be detected also in cr' 
(Fig- EIc)), despite the significant dc contribution: It 
leads to a strong increase at v > 200 MHz, followed by 
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the approach of a plateau close to the upper boundary of 
the investigated frequency range. 




Figure 3: Dielectric-loss spectra of a 5 mmol lysozyme solution in the 
region of the /?- and 5-relaxations at different temperatures. The solid 
Unes are fits using the sum of a Debye function for the /3-relaxation 
and two Cole-Cole functions for the S- and 7-relaxations. Dashed 
lines represent fits with four Debye functions according to Q|. 

Qualitatively similar results as shown in Fig. |2]were 
also obtained for a 5 mmol lysozyme solution. In Fig. [3] 
we show e"(v) within the frequency region of the fi- and 
5-relaxations. Obviously, in this solution with higher 
protein concentration, the 5-relaxation is even more pro- 
nounced with the loss showing clear shoulders close to 
about 200 MHz for the two higher temperatures shown. 

For a meaningful analysis of relaxational processes 
in dielectric spectra, suitable fits, simultaneously per- 
formed for real and imaginary part of the dielectric per- 
mittivity, are essential. In the simplest case, the contri- 
butions of relaxation processes in dielectric spectra can 
be fitted by the Debye equation Ii49il ; 



e*(v) - eoo + 



As 



1 -I- icoT 



(1) 



Ae = fis - fico is the dielectric strength with and 
Soo the limiting values of the real part of the dielectric 
constant for frequencies well below and above the re- 
laxation frequency Vreiax = 1/(2^t), respectively. At 
y - Vreiax a peak shows up in the dielectric loss and 
an inflection point in the frequency dependence of the 
dielectric constant. The dc-conductivity was subtracted 
before the fitting procedure. The Debye theory assumes 
that all dipolar entities relax with the same relaxation 
time T. In reality, however, a disorder-induced distri- 
bution of relaxation times often leads to a considerable 
smearing out of the spectral features llsoi 51 1. An ap- 



Havriliak-Negami formula, which is an empirical exten- 
sion of the Deby e formula by the additional parameters 
gandg Jsllsin: 



e*(v) = £00 + 



Ae 



[1 + (iwT)i-''f 



(2) 



Special cases of this formula are the Cole-Cole formula 
115411 with < < 1 and /? = 1 and the Cole-Davidson 
formula HI 111] with a = and < /3 < 1. While 
the Havriliak-Negami and the Cole-Davidson functions 
are purely empirical, the Cole-Cole distribution of r can 
be approximately derived by the microscopic model of 
a Gaussian distribution of energy barriers, leading to a 
symmetric broadening (compared to the Debye case) of 
the relaxation peak in e" 15711 . 

The y- and y6-relaxations in protein solutions are com- 
monly found to be well describable by Eq. ^ or to 
show at least a behavior very close to monodispersive 
iS, 22, 24, 37] and only few authors apply the Cole- 
Cole function to describe the j6-process 12111 . For the 
y-relaxation, this is reasonable as the corres ponding re - 



propriate phenomenological description is given by the 



laxation of pure water is also of Deb ye ty pe Il48ll58l 5911 
or only slightly broadened ||25, 6^ 61]. Also for the 
y6-relaxation Debye behavior can be expected as any in- 
teraction between the protein molecules is unlikely and 
each molecule "sees" essentially the same environment, 
dominated by the (on the timescale of the j8-relaxation) 
quickly fluctuating water molecules. However, for the 
5-relaxation the situation is far from being clarified, 
especially as the unequivocal detection of its spectral 
shape is hampered by the superposition from the adja- 
cent j6- and y-relaxations (Figs. |2]and[3]). As mentioned 
above, until now it is even not clear if there is only 
a single (5-relaxation or if several relaxation processes 
contribute in this region. In ref. ]0], where lysozyme 
solutions of various concentrations were investigated at 
room temperature, sophisticated arguments favoring the 
use of two Debye functions to describe the (5-relaxation 
at high concentrations were provided. In the present 
work, we fit the experimental data assuming a single 
peak only, which, however, is broadened according to 
the Cole-Cole equation (cf. Eq. (|2]i with /3 = 1). 
When adopting the bound-water explanation of the <5- 
relaxation, a distribution of relaxation times (and thus 
of energy barriers) seems reasonable as there should 
be a variation in the strength of bonding of the water 
molecules of the protein surface. This may depend on 
the polar residue of the macromolecules, to which the 
water molecule is bound, and it may also arise from the 
presence of several hydration shells, the molecules in 
the outer ones being more loosely bound than those in 
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the innermost one |4j. Bound- water relaxations have 
also been previously described by the Cole-Cole func- 



tion 1391 161 
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The Unes in Figs. |2] and [3] are fits with the sum of 
one Debye function for the fi- and two Cole-Cole func- 
tions for the 5- and y-relaxations. Reasonable fits of 
the experimental spectra could be achieved in this way. 
The width parameter as of the (5-relaxation was found to 
vary only weakly around 0. 1 and thus was constrained 
to the range 0.09-0.11 for the final fits. For the y- 
relaxation, the deviation from the Debye case was even 
smaller, the maximum value of ay being a; 0.02. Ac- 
cording to Grant et al. IjtIi . the (5-relaxation of myo- 
globin solutions is bimodal (due to loosely and strongly 
bound water), but the dispersion due to the loosely 
bound water with a peak frequency of a few GHz can be 
incorporated into the y-relaxation by using a Cole-Cole 
function with distribution parameters up to a = 0.07. 
In our case, having even smaller values of a, the use 
of an additional relaxation process {62) seems not jus- 
tified, particularly as an alternative description with a 
Debye function (a - 0) for the y-relaxation does not 
visibly worsen the fits. Nevertheless, for comparison. 
Fig. [3]also shows fits with the sum of four separate De- 
bye functions y, 5\, So) as proposed in ||4|,|37[]. An 
inspection of Fig. [3] by eye does not reveal any signif- 
icant differences in the quality of these fits and, based 
on our data base, we believe it is not possible to make a 
final decision. At least, from the viewpoint of Occam's 
razor, using one relaxation process less for the descrip- 
tion of the data seems preferable. 

The most significant parameter obtained from an 
analysis of relaxational processes is the characteristic 
time of the involved dynamics of the relaxing entities. 
Figure|4]provides the temperature dependence of the re- 
laxation times T of all three detected processes for both 
investigated protein concentrations. All relaxation times 
reveal straight-line behavior in the Arrhenius represen- 
tation of Fig. |4] (dashed lines), indicating thermally ac- 
tivated behavior: 



Toexp 



(3) 



Here, tq is an inverse attempt frequency, often assumed 
to be of the order of a typical phonon frequency and 
denotes the hindering barriers for the relaxational pro- 
cess. However, it should be noted that for all relaxation 
times shown in Fig. |4] deviations from Arrhenius be- 
havior may well be possible when taking into account 
the scatter of the data and the rather small tempera- 
ture range that can be investigated in aqueous solutions, 
which naturally is limited by the freezing and boiling 
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Figure 4: Temperature dependence of the relaxation times obtained 
from the fitting routine. Symbols correspond to t^, Ty, and of the 
3 mmol and 5 mmol lysozyme solution, (a) j8-relaxation; solid lines 
represent Ty, scaled to match rg. (b) (5-relaxation; solid lines represent 
Ty scaled to match r,s. (c) y-relaxation; the solid line corresponds to 
the relaxation times of pure water after |43| . The dashed lines in (a) 
- (c) are linear fits, from which the energy baiTiers were calculated 
according to Eq. {3). 



points of water For pure water, such small deviations 
are well documented 1481 [63ll . In glass forming liquids, 
an often used explanation for non- Arrhenius behavior is 
the cooperativity of the molecular motions 1164 16511 . In 
contrast, for water also a critical power-law of t{T) aris- 
ing from a first-order phase transition was considered 
ll66ll . For comparison, in Fig. Sfc), t{T) of pure water 
from the literature Ii48i1 is included. The absolute val- 
ues of the relaxation times of water are close to those of 
the present y-relaxation times corroborating the assign- 
ment of this relaxation to bulk water molecules. The 
two protein solutions have identical y-relaxation times, 
in accord with the findings in p|. Except for the highest 
temperatures, the relaxation times of pure water seem 
to be slightly lower than those of the protein solutions. 
Slower dynamics of the y-relaxation (i.e. higher values 
of t) than for pur e water was also found in other protein 
solutions i22ll37ll67 l. 
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Whereas the y-relaxation times do not depend on the 
protein concentration, the y6-relaxation times are slightly 
lower for the solution with lower protein concentration 
(Fig. St a)). This finding is consistent with the experi- 
mental results of ref. i22|| for ubiquitin and of |!24] for 
ribonuclease A. It evidences faster reorientational mo- 
tions of the protein molecules in the lower concentrated 
sample, which may be ascribed to weaker hindering of 
this motion by the neighboring molecules. 

Interestingly, irrespective of the absolute values, the 
yS-relaxations of both lysozyme samples show nearly the 
same temperature dependence as the y-relaxation. This 
is visualized in Fig. Bta) by showing Ty{T) (solid lines). 



scaled to match Tp{T) 11681 16911 . Thus, the /^-relaxation 
seems to be strongly coupled to the structural fluctua- 
tions of the solvent (represented by the y-relaxation). 
This is confirmed by the fact that the energy barriers, 
determined from the Arrhenius fits of T/j(T) and Ty{T), 
Eq. (O, are nearly identical (see Fig. |4|. This coupling 
of /3- and y-relaxation simply mirrors the fact that the 
rotation dynamics of a molecule in a solution is essen- 
tially proportional to the solution's viscosity, which is 
expressed by the Debye relation 1271] 



ksT 



(4) 



where 77 is the solvent viscosity and a is the radius of the 
solute molecule. (As the factor l/T in Eq. (|4]l can be 
largely neglected, compared to the exponential or even 
stronger-than-exponential temperature dependence usu- 
ally found for /7(T), in fact Eq. ^ implies oc rj.) 
Just as for pure water jTOil . the viscosity (translational 
motion) and the y-relaxation (reorientational motion of 
the free water molecules) can be assumed to be strongly 
correlated in aqueous solutions and, thus, the finding 
of Tj3 oc Ty (Fig. Eta)) seems reasonable. Equation dU 
in principle allows for the determination of the hydro- 
dynamic radius of the protein molecules, but as they 
are not exactly spherical, this can be a rough approx- 
imation only. As the addition of such small amounts 
of lysozyme does not significantly change the viscos- 
ity of water (which is proofed by the fact that the y- 
relaxation times, which are related to the sample vis- 
cosity, are identical for the solutions and pure water), 
we can use the values for pure water from 117 111 to cal- 
culate the hydrodynamic radius. We arrive at values 
of fl = 1.96 (+0.04) nm for the 5 mmol solution and 
1.91 (+0.02) nm for the 3 mmol sample, largely inde- 
pendent of temperature. This is of same order as the re- 
sults by Bonincontro et al. from dielectric spectroscopy 
(ail.8-l.9nm') 12 ill . Parmaret al. and Chirico et al. from 



light scattering experiments 11721 17311 (1.89+0.025 nm ) 
and Wilkins et al. from Pulse Field Gradient NMR 17411 
(2.05 nm). In the work from Bonincontro et al. 12111 . 
a peak in a(T) was found. Interestingly, in the present 
data for the 5 mmol solution the faint indication of a 
peak is found too, where a varies between 1.93 nm and 
1 .99 nm. In ref. 12111 the peak was ascribed to the tem- 
perature dependence of hydrophobic interactions within 
the protein molecules. Also reversible denaturation ef- 
fects may be considered iItsIi . 

As revealed by Fig. Hb), the (5-relaxation of the 
5 mmol solution is clearly faster than that of the 3 mmol 
solution. Interestingly, in ^ only a single Debye func- 
tion was used for the cJ-relaxation at low concentrations, 
while for the higher ones two were necessary, the sec- 
ond one being located at higher frequencies. Thus, the 
results in may be consistent with the present ones, 
namely the shift of spectral weight to higher frequen- 
cies for higher concentrations, which in our case is di- 
rectly mirrored by the variation of r. A less signifi- 
cant decrease of 5-relaxation times was also observed by 
Oleinikova et al. for Ribonuclease A 12411 . To explain 
this behavior, one could speculate that the structure of 
the water shell around the lysozyme molecule changes 
in dependence of concentration. However, our finding 
that the calculated radius of the lysozyme molecule (in- 
cluding hydration shell) is nearly the same for both con- 
centrations (1.96 vs. 1.91) speaks against such a sce- 
nario. Thus, the microscopic origin of the observed de- 
crease of the ^-relaxation time with increasing concen- 
tration remains unclear. 

The solid lines in Fig. HJb) again show the y- 
relaxation times, scaled to match Tg{l/T). At least for 
the 5 mmol solution, the slope of Tg{l/T) seems to be 
somewhat smaller than for Ty{l/T). This is also mir- 
rored by the energy barriers obtained from the fits with 
Eq. Q (dashed Unes; 0. 17eV for 6 vs. 0.21 eV for y). It 
seems reasonable that the dynamics of the bound water 
molecules is to some extent determined by the interac- 
tions with the polar residues on the protein molecules 
and thus the ^-relaxation is less coupled to the structural 
fluctuations of the solvent, i.e. the y-relaxation. 

The relaxation-time ratio of free and protein-bound 
water dynamics, found in the present work (rs/Ty - 
190 ± 20 for the 3 mmol solution and 130 + 20 for 
5 mmol), is much higher than the factor of about 6-7, 
reported in a recent depolarized light-scattering study 
of lysozyme solutions 1761 . Moreover, the susceptibil- 
ity spectra of that work, extending from 1 GHz well into 
the THz range, reveal significantly faster water dynam- 
ics than commonly detected by dielectric spectroscopy, 
leading to susceptibility peak-frequencies beyond the 
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highest frequencies covered in the present work. As 
noted in ref. ll76ll . these light scattering experiments 
seem to probe mechanisms of different physical origin 
than the molecular reorientations detected by dielectric 
spectroscopy, which points to even more complex water 
dynamics, a fact that certainly warrants further investi- 
gation. 
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Figure 5: Temperature dependence of the inverse relaxation strength 
of the (a) yS-relaxation, (b) i5-relaxation and (c) y-relaxation. Symbols 
correspond to dielectric strengths evaluated by the fitting procedure. 
Lines represent fits according to Eq. (5). The short dashed line in (c) 
coiTesponds to 100/Ae of pure water according to l48h . 

Fig. |5] provides information on the temperature de- 
pendence of the relaxation strengths As of the three ob- 
served relaxation processes. We plot the inverse As, 
which should result in linear behavior for Curie behav- 
ior. 



Ae oc 



(5) 



The latter is expected for dipolar relaxations if the Kirk- 
wood factor, taking into account correlation effects be- 
tween dipoles, is temperature independent l77ll . Indeed 
the relaxation strengths of the 6- and y-relaxations of 
both concentrations can be fitted to Curie-laws (lines in 
Figs. |3b) and (c)). The relaxation strength of the 6- 
relaxation (Fig. |3b)) strongly depends on the concen- 



tration. Assuming bound water as its origin, this is a 
reasonable finding because the number of bound-water 
molecules should increase with the number of proteins. 
Ae^ of the 5 mmol solution is a factor of 1.57 (±0.07) 
higher than for 3 mmol, which is of similar magnitude 
as the 5/3 ratio of the concentrations. This proportion- 
ality is confirmed by the results of ref. 101 (up to a con- 
centration of 1 10 mg/ml, i.e. 7.7 mmol/1), if summing 
up the relaxation strengths of the two <5-relaxations as- 
sumed in this work. 

The absolute values of the strength of the y-relaxation 
(Fig- Etc)) are somewhat smaller than for pure wa- 
ter (see dashed Une in Fig. jS^c) 14811 ) and decrease 
with increasing concentration. This can be primarily as- 
cribed to the trivial substitution effect of water by pro- 
tein molecules in the solution, i.e. the concentration of 
water diminishes with increasing protein concentration. 
As treated in detail in i4J, a further reduction of As is 
caused by a certain amount of water molecules being 
bound to the protein surface, which thus no longer con- 
tribute to the y-relaxation. 

For the y6-relaxation (Fig. |5ja)) the As/^iT) data 
are consistent with Curie behavior at the lower tem- 
peratures, at best. The deviation from Curie behav- 
ior may be ascribed to a temperature variation of the 
dipole moment of the protein molecules as discussed 
below. Alternatively, a temperature-dependent Kirk- 
wood factor may explain the observed deviations. The 
yS-relaxation strength of the 5 mmol solution is by a fac- 
tor of 1.46(+0.05) larger than for 3 mmol, which only 
roughly scales with the expected increase due to the 
larger number density of protein molecules. A similar 
deviation from a purely linear increase of Ae with con- 
centration was also found in 101 for the same concentra- 
tions. It can be explained by a decrease of the effective 
dipole moment jj. of the protein molecules 10, Hi- 

In principle, fi can be calculated from Ae. There are 
a variety of time-honored models enabling such a cal- 
culation, but as a lot of assumptions have to be made 
for these models to be valid, the significance of the ob- 
tained values of fi should not be overrated. Nevertheless 
we used the same approach as in ref. Bl to calculate 
ju by employing the formula predicted by the Onsager- 
Oncley model iHlsoll: 



leoMkBTAs 
Nacsk 



(6) 



Here M is the protein molecular mass, eo is the vacuum 
permittivity, c is the protein concentration in kg/nv' and 
ks and A^a are the Boltzmann and Avogadro constants, 
respectively, denotes the Kirkwood correlation pa- 
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Figure 6: Temperature dependence of the dipole moments fi as cal- 
culated from Eq. j6j for a 3 mmol (triangles) and 5 mmol (circles) 
lysozyme solution. 



rameter, often assumed to be one in diluted protein so- 
lutions |19l|. The obtained results are shown in Fig. |6l 
As expected from the previous paragraph, we obtain 
a lower effective dipole moment for the solution with 
higher protein concentration. In refs. and i24}l this 
was ascribed to antiparallel correlations between differ- 
ent protein molecules, implying gj^ I. The absolute 
values of ju determined from our measurements differ 
from the room-temperature results reported in ref. 101 
(our values are about 1.4 times smaller). This discrep- 
ancy remains unexplained, especially as our findings for 
Ae, used for the calculation of ju, reasonably agree with 
those reported in 10] ■ Takashima et al. have reported 
a room-temperature value (extrapolated to zero concen- 
tration) of 122 Db 14411 . which is of similar magnitude 
as our results. An interesting finding revealed by Fig. 
|6] is the decrease of the dipole moment with tempera- 
ture: While fi{T) is nearly constant for the lower tem- 
peratures for both concentrations, close to room tem- 
perature it starts to decrease with increasing T, which 
becomes especially obvious for the 5 mmol solution, 
for which yU could be determined up to higher tempera- 
tures. This directly mirrors the onset of deviations from 
Curie temperature dependence of As at high tempera- 
tures, documented in Fig. jSja). Proteins are able to 
assume many nearly isoenergetic substates iFtI pTSl pTgll . 
Thus, the observed variation of dipole moment may well 
reflect gradual conformational changes of the molecular 
structure at elevated temperatures. A similar decrease of 
the dipole moment above about 300 K was reported for 
lysozyme solutions (5 mg/ml) of two different pH values 
by Bonincontro et al. 12 lil and attributed to a redistribu- 
tion of microscopic state populations of the protein. In 



that work, a correlation of the temperature dependences 
of fi{T) and the radius a{T) is assumed, while a in the 
present study is nearly temperature independent as men- 
tioned above. 

Fig- 13^) shows the temperature dependence of the 
dc-conductivity cTdc of the two solutions, plotted in a 
way to linearize the Arrhenius behavior predicted for 
ionic conductors, namely 



(7) 



Here, ctq is a prefactor and denotes the hindering 
barrier for the diffusion of the charge carriers. Indeed 
the ionic conductivity closely follows the expected ther- 
mally activated behavior (dashed lines). The fits lead to 
energy barriers of 0.18 and 0.17 eV for 3 and 5 mmol, 
respectively. The conductivity is higher by about a fac- 
tor of two (2.09+0.02) for the 5 mmol solution. An 
increase of the conductivity is reasonable because the 
ions carrying the dc current can be assumed to mainly 
arise from the protein molecules releasing ions when 
dissolved in water. However, it is unclear why the ob- 
served conductivity increase is higher than the concen- 
tration ratio of 5/3. 

In Fig. |7tb) the dc resistivity pdc = l/cdc is shown 
and compared to the scaled y-relaxation times. Espe- 
cially for the 5 mmol sample, clearly different slopes 
of the two curves are revealed; indications for similar 
behavior is also found for 3 mmol. This finding im- 
plies a breakdown of the Debye-Stokes-Einstein rela- 
tion, which can be expressed as p^c ^ Ty and is an often 
discussed phenomenon in supercooled liquids 18^ 81]. 
Interestingly this decoupling of charge transport and y- 
relaxation closely resembles the one between the S- and 
y-relaxation documented in Fig. ISb). Moreover, the 
energy barriers for dc transport (Fig. 121^)) and for the 6- 
relaxation (Fig. 3(b)) are nearly identical. This surpris- 
ing finding is difficult to rationalize and may be acciden- 
tal. However, one should be aware that the conductivity 
is proportional to both, the mobiUty and the number of 
charge-carrier. If one assumes that the temperature de- 
pendence of (TdciT) in the present case is dominated by 
the number of ions released from the protein molecules 
rather than the mobility, a close connection of both pro- 
cesses seems possible: For the 5-relaxation, i.e. the 
reorientational motions of bound water molecules, the 
bonds to the protein molecules have to be (temporar- 
ily) broken, which may be determined by similar energy 
barriers as those necessary for the release of ions into 
the solution. Further work is necessary to corroborate 
this speculation, e.g., by investigating the relation of 
conductivity and i5-relaxation for other protein species. 
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Figure 7: (a) Arrhenius type presentation of (a) tlie dc-conductivity 
CTjc of a 3 mmol (triangles) and a 5 mmol (circles) lysozyme solution 
and (b) the dc-resistivity pjc = l/iTdc- The dashed lines are linear fits 
to the data, corresponding to Eq. Q. The solid fines represent scaled 
values of Ty to test for the Debye-Stokes-Einstein relation, pj^ Ty. 



4. Conclusions 

In the present work, we have provided a thorough 
characterization of two protein solutions with differ- 
ent concentrations using high-frequency dielectric spec- 
troscopy from 1 MHz to 40 GHz. A variety of informa- 
tion on the j0- and y-relaxations arising from the protein 
tumbling and reorientation of the free water molecules 
has been collected. Most importantly we have detected 
a well-pronounced (5-dispersion, attributed to bound wa- 
ter dynamics, and have obtained detailed information on 
its temperature dependence. Using a Debye function for 
the y6 process, a Cole-Cole function for the y-relaxation, 
and a single Cole-Cole function for the description of 
the (5-dispersion, the complete broadband spectra can be 
well fitted. 

Temperature-dependent data on the relaxation time 
and strength have been obtained for aU three main dis- 
persion regions of lysozyme solutions enabling the de- 
termination of hindering barriers for the relaxational 
processes and for the diffusion of ionic charge carriers. 
Obviously all energy barriers in these protein solutions 
are of similar order of magnitude, varying between 0. 17 
and 0.22 eV. While we find the expected strong correla- 
tion of the jS- and y-relaxation, the (5-relaxation seems to 
be less strongly influenced by the fluctuations of the sol- 
vent and, instead, is governed by interactions with the 



protein molecules. We have found a significant concen- 
tration dependence of the (5-relaxation dynamics, whose 
origin is unclear until now. From our results we have 
deduced the hydrodynamic radius and the temperature 
dependence of the dipole moment. A noticeable result is 
the decrease of the latter with increasing temperatures, 
which also leads to deviations of the relaxation strength 
of the j8-relaxation from Curie behavior. We attribute 
this finding to gradual conformational changes of the 
protein structure. 

The analysis of the temperature-dependent dc con- 
ductivity and its comparison with the y-relaxation time 
reveals a breakdown of the Debye-Stokes-Einstein rela- 
tion, i.e. the ionic charge transport is governed by dif- 
ferent energy barriers than the motions of the solution 
molecules. Interestingly, we find that instead the charge 
transport and the 5-relaxation, i.e. the reorientation of 
bound water molecules, are determined by identical en- 
ergy barriers, an unexpected and so far unexplained be- 
havior 

Overall, our high-frequency dielectric measurements 
demonstrate the rich dynamics of protein solutions, 
which shows many properties about whose microscopic 
origins currently only speculations are possible. It is 
clear that further work is needed, especially covering 
a broad frequency range and involving temperature- 
dependent measurements. 
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